Radiation-induced magnetotransport in high-mobility two-dimensional systems: 

Role of electron heating 
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Effects of microwave radiation on magnetoresistance are analyzed in a balance-equation sclieme 
that covers regimes of inter- and intra-Landau level processes and takes account of piioton-asissted 
electron transitions as well as radiation-induced change of the electron distribution for high mo- 
bility two-dimensional systems. Short-range scatterings due to background impurities and defects 
are shown to be the dominant direct contributors to the photoresistance oscillations. The electron 
temperature characterizing the system heating due to irradiation, is derived by balancing the en- 
ergy absorption from the radiation field and the energy dissipation to the lattice through realistic 
electron-phonon couplings, exhibiting resonant oscillation. Microwave modulations of Shubnikov de 
Haas oscillation amplitude are produced together with microwave-induced resistance oscillations, in 
agreement with experimental findings. In addition, the suppression of the magnetoresistance caused 
by low- frequency radiation in the higher magnetic field side is also demonstrated. 

PACS numbers: 73.50.Jt, 73.50.Mx, 78.67.-n, 78.70.-g 



I. INTRODUCTION 



The discovery of microwave induced magnetor- 
ersistance oscillations (MIMOs) and zero-resistance 
states (ZRS) in high-mobility two-dimensional 
(2D) electron gas (EG)^'^'^''' has stimulated 
tremendous experimental^iSiLSiSiiSiiiiiSiiiii and 
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interest in radiation related magneto-transport in 2D 
electron systems. Since theoretically it has been shown 
that the ZRS can be the result of the instability in- 
duced by absolute negative resistivityji^ the majority 
of microscopic models focus mainly on MIMOs in 
spatially uniform cases and identify the region where 
an negative dissipative magnetoresistance develops as 
that of measured zero resistance. Most of previous 
investigations concentrated on the range of low magnetic 
fields uJc/uj < 1 {ujc stands for the cyclotron frequency) 
subject to a radiation of frequency uj/2tt < 100 GHz, 
where MIMOs show up strongly and Shubnikov-de 
Haas oscillations (SdHOs) are rarely appreciable. In 
spite of the fact that both MIMOs and SdHOs are 
magnetoresistance related phenomena appearing in 
overlapping field regimes, little attention was paid to 
the influence of a microwave radiation on SdHO until 
a recent experimental finding at higher frequencyiii 
Further observations clearly show that the amplitudes 
of SdHOs are strongly affected by microwave radiations 
of different frequency in both low {u>c/lo < 1) and 
high {iVc/io > 1) magnetic field ranges liSiiiii Kovalev 
et alm^ observed a suppression of the SdHOs around 
cyclotron resonance Wc ~ w induced by a radiation of 
285 GHz. Du et aZjiS found strong modulations of SdHO 
in an ultra-clean 2D sample subjected to microwaves of 
146 GHz, clearly showing, in addition to the first node 
at uj/ujc = 1, higher order nodes around uj/uJc = 2 and 3. 
Mani— reported strong modulation in the amplitude of 
SdHOs accompanying MIMOs and zero-resistance states 



excited by a 163.5-GHz radiation and large dropoff of 
the dissipative resistivity below its dark value at high 
{ujc/uj > 1) field side when subjected to low- frequency 
radiation. Very recently, Dorozhkin et alm^ reported 
both the strong suppression of the magnetoresistance 
caused by radiation below 30 GHz and an interesting 
modulation of SdHOs in the range lUc/'-^ > 1- They 
found that SdHOs are generally strongly damped by the 
radiation but there is a narrow magnetic field range in 
between allowed ranges of inter- and intra-Landau level 
transitions, where the amplitude of SdHO is insensitive 
to the microwave irradiation. These observations pro- 
vide a more complicated and appealing picture of the 
microwave-related transport phenomena, which must be 
accounted for in any theoretical model for MIMOs. 

We propose that these SdHO modulations come from 
the electron heating induced by the microwave radiation. 
Under the illumination of microwave the electron system, 
which continuously absorbs energy from the radiation 
field, would certainly be heated. Unfortunately, the elec- 
tron heating has so far been ignored in most of the theo- 
retical treatments. The electron-acoustic phonon interac- 
tion was previously considered to contribute to Landau- 
level broadeningSi or to act as a dampingSi for the orbit 
movement, providing a mechanism for the suppression 
of MIMOs when the lattice temperature increases. Be- 
sides the inelastic electron-phonon scattering also plays 
another important role to dissipate energy from the elec- 
tron system to the lattice. The energy absorption rate 
is indeed small in high-mobility electron systems at low 
temperature as in the experiments. This, however, does 
not imply a negligible electron heating, since the elec- 
tron energy-dissipation rate is also small because of weak 
electron-phonon scattering at temperature T < 1 K. To 
deal with SdHO, which is very sensitive to the smearing 
of the electron distribution, one has to carefully calcu- 
late the electron heating due to microwave irradiation in 
a uniform model. 

On the other hand, microwave irradiation heats the 
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electrons and thus greatly strengthens the thermalizing 
trend of the system by enhancing the electron-electron 
scattering rate at this low temperature regime. This en- 
ables us to describe these high-mobility 2D electron sys- 
tems with a quasi-equilibrium distribution in a moving 
reference frame. 

In this paper we pursue a theoretical investigation on 
MIMOs and SdHOs taking account of the electron heat- 
ing under microwave irradiation. We generalize the bal- 
ance equation approach to radiation-induced magneto- 
transport in high mobility two-dimensional electron sys- 
tems. By carefully calculating the electron heating based 
on the balance of the energy absorption from the ra- 
diation field and the energy dissipation to the lattice 
through electron-phonon interactions in a typical GaAs- 
based heterosystem and taking into account the electro- 
dynamic effect, we are able not only to reproduce the 
interesting phenomena of MIMOs in quantitative agree- 
ment with experiments in amplitudes, phases and radi- 
ation dependence of the oscillation, but also to obtain 
SdHO modulations observed in the experiments. 



II. FORMULATION 
A. Force- and energy-balance equations 

This paper is concerned with the magnetotransport in 
a microscopically homogeneous 2D system, and refers the 
measured zero resistance to the macroscopic consequence 
of the instability due to the occurrence of negative dissi- 
pative resistivity. 

We consider A^e electrons in a unit area of an infinite 
quasi-2D system in the x-y plane with a confining poten- 
tial V{z) in the z direction. These electrons, in addition 
to interacting with each other, are scattered by random 
impurities and/or disorders and by phonons in the lat- 
tice. Within the magnetic field range relevant to MIMO 
phenomenon, the experiments exclude the onset of the 
quantum Hall effect, thus allowing us to assume that the 
2D electrons are in extended states. 

To include possible elliptically polarized microwave il- 
lumination we assume that a dc electric field Eq and a 
high-frequency (HF) ac field of angular frequency lu, 

E(t) EE E, sin(wt) + Eccos(wi), (I) 

are applied inside the 2D system in the x-y plane, to- 
gether with a magnetic field B = (0, 0, B) along the z 
direction. The spatial homogeneity of the fields and the 
parabolic band structure allows to describe the trans- 
port of this system in terms of its center-of-mass (cm.) 
motion and the relative motion, i.e. the motion of elec- 
trons in the reference frame moving with the c.m^^'^^^ 
The center-of-mass momentum and coordinate of the 2D 
electron system are defined as P = J2jPj\\ ^^"^ ^ = 
^e^J2j^j\\ with PjIi = {pjx,Pjy) and = {xj,yj) be- 
ing the momentum and coordinate of the jth electron in 



the 2D plane, respectively, and the relative electron mo- 
mentum and coordinate are defined as p^|| = p^n — P/A'o 
and r^ii = rj|| — R, respectively. In terms of these vari- 
ables, the Hamiltonian of the system, H, can be written 
as the sum of a center-of-mass part Hem and a relative 
electron part Her (A(r) is the vector potential of the B 
field), 

i?cm - — ^(P-AeeA(R))2-iV„e(Eo + E(<)).I(?) 

^^•-E[i(p;„-A(r;,))V|j--,y(.,)] 
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+ ^yc(r-|| -r^l|,z^,2;j), (3) 

together with electron-impurity and electron-phonon in- 
teractions 

Tlei = E "(qi|,^a)e'''"-^''+'"^" \ (4) 
Fep = ^M(q,|,g.)(feq-f6t)e^<i,r(i^+r^„). (5) 

Here m and are, respectively, the electron effective 
mass parallel and perpendicular to the 2D plane, and 
Vc stands for the electron-electron Coulomb interaction; 
M(q||,Z(j) is the potential of the ath impurity locating 
at {ra\\,Za); ^q(&q) are the creation (annihilation) oper- 
ators of the bulk phonon with wavevector q = (q||,(7z) 
and M(q||,(72) is the matrix element of the electron- 
phonon interaction in the 3D plane-wave representation. 
Note that the uniform electric field (dc and ac) appears 
only in i?cm, and that Her is just the Hamiltonian of 
a quasi-2D system subjected to a magnetic field with- 
out an electric field. The coupling between the center- 
of-mass and the relative electrons appears only in the 
exponential factor exp(iq|| • R) inside the 2D momem- 
tum qy summation in Hd and Hcpi^ The balance equa- 
tion treatment starts with the Heisenberg operator equa- 
tion for the rate of change of the center-of-mass veloc- 
ity V = -i[V,H] + dW/dt with V = -i[R,i?], and 
that for the rate of change of the relative electron energy 
Her = —'AMb-ctH]. Then we proceed with the determina- 
tion of their statistical averages. 

As proposed in Ref.l35l the cm. coordinate opera- 
tor R and velocity operator V can be treated classi- 
cally, i.e. as the time-dependent expectation values of 
cm. coordinate and velocity, R(i) and V(t), such that 
R(t) - R(i') = Y{s)ds. We are concerned with the 
steady transport state under an irradiation of single fre- 
quency and focus on the photon-induced dc resistivity 
and the energy absorption of the HF field. These quanti- 
ties are directly related to the time-averaged and/or base- 
frequency oscillating components of the cm. velocity. 
Although higher harmonics of the current may affect the 
dc and lower harmonic terms of the drift velocity through 
entering the damping force and energy exchange rates 
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in the resulting equations, in an ordinary semiconductor 
the power of even the third harmonic current is rather 
weak as compared to the fundamental current. For the 
HF field intensity in the MIMO experiments, the effect 
of higher harmonic current is safely negligible. Hence, it 
suffices to assume that the cm. velocity, i.e. the electron 
drift velocity, consists of a dc part vq and a stationary 
time-dependent part v{t) of the form 



V(t) = vo — vi cos{Lut) — V2 sin(wi). 



(6) 



With this, the exponential factor in the operator equa- 
tions can be expanded in terms of Bessel functions J„(a;), 

nt °° 



n— — oo 



771^0 n= — co 



Here the argument in the Bessel functions 



1 



^^-[(q||.Vi)^ + (q||-V2)2]% 



(7) 



and tan.^ = (q|| • V2)/(q|| • Vi). 

Under the influence of a modest-strength HF electric 
field the electron system is far from equilibrium. How- 
ever, the distribution function of relative electrons, which 
experience no electric field directly, may be close to an 
quasi-equilibrium type distribution function. For the 
experimental GaAs-based ultra-clean 2D electron sys- 
tems having carrier mobility of the order of 2000m2/Vs, 
the elastic momentum scattering rate is around ~ 
10 mK. In these systems, the thermalization time rth (i.e. 
the time for system to return to its internal equilibrating 
state when it is deviated from) , estimated conservatively 
using electron-electron (e-e) interaction related inelastic 
scattering time Tee calculated with an equilibrium distri- 
bution function at temperature T = IK, is also around 
Tj^^ ^ Tj^^ ^ 10 niK. The illumination of microwave 
certainly heats the electrons. Even an electron heat- 
ing comparable to a couple of degrees temperature rise 
would greatly enhance T~e^ , such that the thermalization 
time Tth would become much shorter than the momen- 
tum relaxation time Tm under microwave irradiationi^ 
The relative electron systems subject to a modest radia- 
tion would rapidly thermalize and can thus be described 
reasonably by a Fermi-type distribution function at an 
average electron temperature Te in the reference frame 
moving with the center-of-mass. This allows us to carry 
out the statistical average of the operator equations for 
the rates of changes of the cm. velocity V and relative 
electron energy Her to the leading order in Hd and H^p 
with succinct forms. 

For the determination of unknown parameter Vg, Vi, 
V2, and Te, it suffices to know the damping force up to the 
base frequency oscillating term F{t) = Fq -I- Fg sin{ujt) + 
FcCos{ujt), and the energy-related quantities up to the 



time-average terms. We finally obtain the force and en- 
ergy balance equations, 



A^ceEo + iVce(vo x B) + Fq = 0, 



(8) 



eE, 



vi 



V2 



moj NciriLU muj 
eEc Fr e 



rriLO Ncinuj muj 



(V2 X B), (9) 

(vi X B), (10) 
(11) 



Here 



Fo = ^|[/(q||)f q||^'(0n2(q||,c.o-nt^) 

qy 71— — OQ 

OO 

+ ^|Af(q)|' J2 q||J,^(OA2(q,c.o + f^q-nc^)(12) 

q n— — oo 

is the time-averaged damping force, 

oo 

qy n— — oc 

oo 

+ ^|M(q)|' nu; J^{0A2{CL, coo + n^~nu;) {13) 



is the time-averaged rate of the electron energy absorp- 
tion from the HF field, and 



M/ = ^|M(q)|2 Y n^J^{OA2{ci,u;o + n^-nu;) 



n— — OQ 



(14) 

is the time-averaged rate of the electron energy dissi- 
pation to the lattice due to electron-phonon scatterings. 
The oscillating frictional force amplitudes F^ = F22 — Fn 
and Fc = F21 + F12 are given by (/i = 1, 2) 

00 

Fi^ = -^q||r7^|t/(q)l)P ^ [J^O]' ^i(^\\,^o - nco) 

qy n— — 00 

oc 

-^q||r;^|M(q)p ^ Ai{ci,UJo + - nu;),{15) 



n— — OC 



F2,=Eqil7l^(q||)l' ^ 2nJ^(On2(q||,c.o-r.a;) 

qy 11 — ~oo 

00 

+ ^q||^|Af(q)p Y 2nj2(OA2(q,c.o + f^q-nc^). (16) 

q n— — oc 

In these expressions, 77^ = qy • v^/w^; luq = 
q|| • Vq; ?7(q||) and M(q) are effective impurity and 
phonon scattering potentials (including effects of the 
spatial distribution of impurities and the form fac- 
tor of quasi-2D electrons)i^ n2(q||,ri) and A2(q, fi) = 
2H2(q||,rJ)[n(r2q/r)-n(r2/re)] (withn(a;) = l/(e^-l)) 
are the imaginary parts of the electron density correlation 
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function and electron-phonon correlation function in the 
presence of the magnetic field. ni(q||,r2) and Ai(q, fi) 
are the real parts of these two correlation functions. 

Effects of a microwave radiation on electron transport 
first come from the HF field induced cm. motion (elec- 
tron drift motion) and the related change of the electron 
distribution. In addition to this, the HF field also enters 
via the argument ^ of the Bessel functions in Fq, F^^, W 
and Sp. Compared with that without a HF field, we see 
that in an electron gas having impurity and/or phonon 
scatterings (otherwise homogeneous), a HF field of fre- 
quency Lo opens additional channels for electron transi- 
tion: an electron in a state can absorb or emit one or sev- 
eral photons of frequency uj and scattered to a different 
state with the help of impurities and/or phonons. The 
sum over |n| > 1 represents contributions of real single 
and multiple photon participating processes. The role of 
these processes is two folds. On the one hand, they con- 
tribute additional damping force to the moving electrons, 
giving rise directly to photoresistance, and at the same 
time, transfer energy from the HF field to the electron 
system, resulting in electron heating, i.e. another change 
(smearing) in the electron distribution.^^ Furthermore, 
the radiation field, showing up in the term with Jo{^) in 
Fo, F^i/ and W, gives rise to another effective change of 
damping forces and energy-loss rate, without emission or 
absorption of real photons. This virtual photon process 
also contributes to photoresistancei^ All these effects are 
carried by parameters Vq, Vi, V2 and T^. Eqs. l|5|l- Hll|) 
form a closed set of equations for the determination of 
these parameters when Eg, Ec and Eg are given in a 2D 
system subjected to a magnetic field B at temperature 
T. 



B. Longitudinal and transverse resistivities 

The nonlinear resistivity in the presence of a high- 
frequency field is easily obtained from Eq. ijHJi. Taking 
Vq to be in the x direction, vq = (uox,0, 0), we immedi- 
ately get the transverse and longitudinal resistivities, 



Rxx 



En 



NcevQx NcB 

Eqx Fq 



(17) 
(18) 



The linear magnetoresistivity is the weak dc current limit 
{vox 0): 



n— — c« 



n— — 00 



(19) 



Note that although according to Eas. (|12|l . (|18|) and 
(|19|l . the longitudinal magnetoresistivity Rxx can be for- 



mally written as the sum of contributions from various in- 
dividual scattering mechanisms, all the scattering mech- 
anisms have to be taken into account simultaneously in 
solving the momentum- and energy-balance equations 
@, ((TUI) and |[TT|) for vi, V2 and To, which enter the 
Bessel functions and other parts in the expression of Rxx- 



C. Landau- level broadening 

In the present model the effects of interparticle 
Coulomb screening are included in the electron com- 
plex density correlation function H(q|| , J7) = Hi(q|| , ft) + 
«H2(q||,r2), which, in the random phase approximation, 
can be expressed as 



H(q||,17) 



no(q||,ri) 



where 



e(q||,17) = l-y(gt|)Ho(q||,r!) 
is the complex dynamical dielectric function. 



^(911 



2eoKq 



(20) 



(21) 



(22) 



is the effective Coulomb potential with k the di- 
electric constant of the material and H{q\\) is 
a 2D wavefunction-related overlapping integration;^ 
Ho(q||, ri) — Hoi(q||, ri) + iHo2(q||, fi) is the complex den- 
sity correlation function of the independent electron sys- 
tem in the presence of the magnetic field. With this dy- 
namically screened density correlation function the col- 
lective plasma modes of the 2DES are incorporated. Dis- 
regard these collective modes one can just use a static 
screening e(q||,0) instead. 

The Ho2(q||, ^) function of a 2D system in a magnetic 
field can be written in terms of Landau representations^ 

Ho2(q|| .^) = 7^Y1 al'?'/2)H2(n, n', 17), (23) 



H2(n,n',r2) 



de[f[e)- f{e + n)] 



X ImG„(£ + f7)ImG„'(£), (24) 



where Ib ~ ■\/l/|ei?| is the magnetic length. 



71, 



{Y)^v\[{n + m-'Y'e-^[Li{Y)] 



(25) 



with L^n(X) the associate Laguerre polynomial, f{e) = 
{exp[(e — ii)/Tc] -\- 1}~^ the Fermi distribution function, 
and ImGn (e) is the imaginary part of the electron Green's 
function, or the density of states (DOS), of the Landau 
level n. The real part function Hoi(q||,r2) and corre- 
sponding Aoi(q||,f2) function can be derived from their 
imaginary parts via the Kramers-Kronig relation. 
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In principle, to obtain the Green's function ImG„(e), a 
self-consistent calculation has to be carried out from the 
Dyson equation for the self-energy with all the impurity, 
phonon and e-e scatterings included. The resultant G„ 
is generally a complicated function of the magnetic field, 
temperature, and Landau-level index n, also dependent 
on the different kinds of scatterings. Such a calculation 
is beyond the scope of the present study. In this paper 
we model the DOS function with a Gaussian-type form 
(£„ is the energy of the n-th Landau level) i^2i^ 



27r 



ImG„(e) = j:^cxp 



Y2 



with a broadening width given by 

/^eujra\ 1/2 

r= — -) , 



(26) 



(27) 



where /xo is the linear mobility in the absence of the mag- 
netic field and a is a semiempirical parameter to take into 
account the difference of the transport scattering time 
Tm determining the mobility hq, from the single parti- 
cle lifetime Tg related to Landau level broadening. The 
latter depends on elastic scatterings of different types 
and their relative strengths, as well as contributions of 
electron-phonon and electron-electron scatterings, a will 
be served as the only adjustable parameter in the present 
investigation. Unlike the semielliptic function, which 
can model only separated Landau-level Gaussian- 
type broadening function can reasonably cover both the 
separated-level and overlapping-level regimes. 



D. Effect of radiative decay 

The HF electric field E(t) appearing in Eqs. (8) and 
(9) is the total (external and induced) field really acting 
on the 2D electrons. Experiments are always performed 
under the condition of giving external radiation. In this 
paper we assume that the electromagnetic wave is inci- 
dent perpendicularly (along z-axis) upon 2DEG from the 
vacuum with the incident electric field of 



Ei(i) = Eiis siii{ujt) + Eic cos{Lut) 



(28) 



at plane z = 0. The relation between E(i) and Ei(t) 
is easily obtained by solving the Maxwell equations con- 
necting both sides of the 2DEG which is carrying a sheet 
current density Ncev{t). If the 2DEG locates under the 
surface plane at z = of a thick (treated as semi-infinite) 
semiconductor substrate having a refraction index n^, we 
have^ 



E(t) = 



iVcev(i) 2 



-Ei(t). 



(29) 



If the 2DEG is contained in a thin sample suspended in 
vacuum at the plane z = 0, then 



2eoc 



Ei(t). 



(30) 



In the numerical calculation of this paper we consider 
the latter case and use Eq. (|3U|1 for the total selfconsis- 
tent field E(i) in Eqs. © and ifTIUl . This electrodynamic 
effect 4^*^ recently refered as radiative decay,^'' gives rise 
to an additional damping in the 2DEG response to a 
given incident HP field. The induced damping turns out 
to be much stronger than the intrinsic damping due to 
scattering-related forces Fg and Fc for the experimen- 
tal high-mobility systems at low temperatures. For al- 
most all the cases pertinent to MIMO experiments we 
can neglect the forces F^ and Fc completely in solving 
vi = (vixjViy) and V2 = {v2x, V2y) from Eqs.Q and l(Tn|) 
for given incident fields Ei<; and Ejc, and obtain explicitly 



Vlx = (aXsx + bXsy)/^ 
Vly = {aXsy - bXsx)/^ 

V2x = {-axcx-bxcy)/^ 

V2y = {-aXcy + bXcx)/^ 

with A ={1-51+ 7')' + (27c.^c.)^ and 



(31) 



Xsx 
Xsy 
Xcx 
Xcy 



+ 5^Vcx + IloV, 
+ S^^sy - lu. 
l^cy - 5ujVsx + lujVs' 



cy 



(32) 



Here 



LOc/uj and = ^/lu with 

2meoc 



(34) 



With these Vi and V2, the argument ^ entering the Bessel 
functions is obtained. All the transport quantities, such 
as Sp, W and Rxxi can be calculated directly with the 
electron temperature To determined from the energy bal- 
ance equation ((TT)l . 



III. NUMERICAL RESULTS FOR GAAS-BASED 
SYSTEMS 

As in the experiments, we focus our attention on high 
mobility 2DEGs formed by GaAs/AlGaAs heterojunc- 
tions. For these systems at temperature T < IK, the 
dominant contributions to the energy absorption Sp and 
photoresistivity Rxx — RxxiO) come from the impurity- 
assisted photon-absorption and emission process. At 
different magnetic field strength, this process is asso- 
ciated with electron transitions between either inter- 
Landau level states or intra-Landau-level states. Ac- 
cording to H26|l . the width of each Landau level is about 
2r. The condition for inter-Landau level transition with 
impurity-assisted single-photon process^"' is ui > ujc — 2r, 
or ujc/uj < flintcr = (/? + \/ [3'^ + 4)^/4; and that for 
impurity-assisted intra-Landau level transition is u; < 2r, 
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or LOcI'-^ > ciintra = here f3 = (32ea / TrmfiQUj)^ . 

However, since the DOS of each Landau level is assumed 
to be Gaussian rather than a clear cutoff function and 
the multi-photon processes also play roles, the transition 
boundaries between different regimes may be somewhat 
smeared. 

As indicated by experimentSj*- although long range 
scattering due to remote donors always exists in the 
2D heterostructures, in ultra-clean GaAs-based 2D sam- 
ples having mobility of order of 10"^ m^/Vs, the remote 
donor scattering is responsible for merely 10% or less 
of the total momentum scattering rate. The dominant 
contribution to the momentum scattering rate comes 
from short-range scatterers such as residual impurities 
or defects in the background. Furthermore, even with 
the same momentum scattering rate the remote impu- 
rity scattering is much less efficient in contributing to 
microwave-induced magnetoresistance oscillations than 
short-ranged background impurities or defects,'*^ There- 
fore, in the numerical calculations in this paper we as- 
sume that the elastic scatterings are due to short-range 
impurities randomly distributed throughout the GaAs 
region. The impurity densities are determined by the 
requirement that electron total linear mobility at zero 
magnetic field equals the giving value at lattice temper- 
ature T. Possibly, long-range remote donnor scattering 
may give rise to important contribution to the Landau- 
level broadening. This effect, together with the role of 
electron-phonon and electron-electron scatterings, is in- 
cluded in the semiempirical parameter a in the expression 

In order to obtain the energy dissipation rate from 
the electron system to the lattice, W, we take into ac- 
count scatterings from bulk longitudinal acoustic (LA) 
and transverse acoustic (TA) phonons (via the defor- 
mation potential and piezoelectric couplings), as well as 
from longitudinal optical (LO) phonons (via the Frohlich 
coupling) in the GaAs-based system. The relevant ma- 
trix elements are well known^ The material and cou- 
pling parameters for the system are taken to be widely 
accepted values in bulk GaAs: electron effective mass 
m = 0.068 rrie (^tIg is the free electron mass), transverse 
sound speed Vst = 2.48 x lO^m/s, longitudinal sound 
speed Vsi = 5.29xl0'^m/s, acoustic deformation potential 
S = 8.5 eV, piezoelectric constant eu — 1.41 x 10^ V/m, 
dielectric constant k — 12.9, material mass density d = 
5.31 g/cur^. 



100 GHz 



Figure 1 shows the calculated energy absorption rate 
Sp, the electron temperature To and the longitudinal 
magnetoresistivity Rxx as functions of Wc/w for a 2D sys- 
tem having an electron density of — 3.0 x 10^^ m~^, 
a linear mobility of fig = 2000ni-^/Vs and a broad- 
ening parameter of a = 10, subject to linearly x- 
direction polarized incident microwave radiations of fre- 
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FIG. 1: The magnetoresistivity -R^a;, electron temperature 
Te and energy absorption rate of a GaAs-based 2DEG 
with /io = 2000 mVVs and a = 10, subjected to 100 GHz 
linearly a;-polarized incident HP fields E\a sin(ajt) of four dif- 
ferent strengths. The lattice temperature \&T =\K. 



100 GHz T=1K 

P.= 210 W/m' 



w/o radiation 
linear x-pol 
linear y-pol 
left circular 
right circular 




FIG. 2: The magnetoresistivity K:cx versus w/tJc for the same 
system as described in Fig. 1, subject to 100 GHz HF radi- 
ations of fixed incident power of Pi = 210 W/m^ but four 
different polarizations. 



quency u) jli: — 100 GHz having four different amplitudes 
E\s = 2.2,3,4 and 5 V/cm at a lattice temperature of 
T = 1 K. The energy absorption rate 5p exhibits a broad 
main peak at cyclotron resonance LUc/to = 1 and sec- 
ondary peaks at harmonics lUc/^ — 1/2,1/3,1/4. The 
electron heating has similar feature: To exhibits peaks 
around cUc/uj = 1,1/2,1/3,1/4. For this GaAs system 
P = 0.65, aintor = 1-6 and ai„tra — 4.7. We can see 
that, at lower magnetic fields, especially oJc/uj < 1.4, the 
system absorbs enough energy from the radiation field 
via inter-Landau level transitions and Tc is significantly 
higher than T, with the maximum as high as 21 K around 
LUc/i^ — 1. With increasing strength of the magnetic 
field the inter-Landau level transition weakens (impurity- 
assisted single-photon process is mainly allowed when 
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ijJcl^ < flintcr = 1-6) and the absorbed energy decreases 
rapidly. Within the range 2 < Wc/w < 4 before intra- 
Landau level transitions can take place, Sp is two or- 
ders of magnitude smaller than that in the low magnetic 
field range. Correspondingly the electron temperature 
To is only slightly higher than the lattice temperature T. 
The magnetoresistivity R^x showing in the upper part of 
Fig. 1, exhibits interesting features. MIMOs (with fixed 
points rather than extrema at uJc/(^ = 1, 1/2, 1/3, f /4) 
clearly appear at lower magnetic fields, which are in- 
sensitive to the electron heating even at of order of 
20 K. SdHOs appearing in the higher magnetic field side, 
however, are damped due to the rise of the electron 
temperature > 1 K as compared with that without 
radiation. With an increase in the microwave ampli- 
tude from Eis = 2.2V/cm to 5V/cm, MIMOs become 
much stronger and SdHOs are further damped. But 
the radiation-induced SdHO damping is always relatively 
smaller within 2.4 < uJc/^^ < 4 between allowed ranges of 
inter- and intra-Landau level transitions. 

It is worth noting that the predicted MIMOs here ex- 
hibit much improved agreement with experiments over 
previous theoretical models. The maxima of R^x oscilla- 
tion locate at lu/uJc = and minima at lo/lOc = j+S+, 
with(5± ~ 0.23-0.25 for j = 2, 3, 4... and J± ~ 0.16-0.18 
for j = I (see Fig. 2). These phase details, as well as the 
absolute (rather than reduced) magnitudes of the oscil- 
lation amplitudes and the required incident microwave 
strengths to induce oscillations are in good quantitative 
agreement with experiments^s^i^i^ 

The MIMOs depend on the polarization of the inci- 
dent microwave field in respect to the dc field Eq. Physi- 
cally this is clear in the present model since it is through 
the cm. motion that a HF field affects the photore- 
sistivity of the 2D electron system. Under the influ- 
ence of a magnetic field perpendicular to the plane, the 
cm. performs a cyclic motion of frequency cOc in the 2D 
plane. A perpendicularly incident circularly-polarized 
microwave would accelerate or decelerate this cyclic mo- 
tion depending on the HF electric field circling with or 
against it. Thus, at fix incident power, a left-polarized 
microwave would yield much stronger effect on the Rxx 
oscillation than a right-polarized one and this effect is ap- 
parently strongest in the vicinity of cylcotron resonance 
Wc/w = 1. The difference between the x-direction lin- 
early polarized wave and the y-direction linearly polar- 
ized wave, however, comes mainly from the the differ- 
ent angle of radiation-induced cm. motion with respect 
to the dc current, and thus not so sensitive to that of 
the uJc/^ range. In Fig. 2 we plot the calculated Rxx 
versus uj/lJc for the same system as described in Fig. 1, 
subject to a 100 GHz microwave radiation having a fixed 
incident power of Pi = 210W/m2 (equivalent to an inci- 
dent amplitude Eis — 4V/cm of linear polarization) but 
four different polarizations: linear x-polarizaton, linear 
y-polarization, left circular polarization and right circu- 
lar polarization. Their difference is clearly seen. 



E.^ = 0.8 V/cm 
E^^ = 1 .2 V/cm 
e!J = 2.0 V/cm 
bI = 3.5 V/cm 




CO /o) 



FIG. 3: The magnetoresistivity Rxx, electron temperature 
Tc and energy absorption rate Sp of a GaAs-based 2DEG 
with fio = 2500m'^/Vs and a — 12.5, subjected to 50 GHz 
linearly x-polarized incident HF fields Eia sin(a;t) of four dif- 
ferent strengths. The lattice temperature is T = 1 A". 




FIG. 4: The magnetoresistivity Rxx versus uj/lJc for the same 
system as described in Fig. 3, subject to 50GHz linearly x- 
polarized incident HF fields of four different strengths. 



B. 50 GHz and lower frequency 

Figure 3 shows the energy absorption rate S'p, the elec- 
tron temperature Te and the longitudinal magnetoresis- 
tivity Rxx as functions of Wc/w for a 2D system having 
an electron density of A^c — 3.0 x 10^^ m^^, a linear mo- 
bility of fiQ = 2500m^/Vs and a broadening parameter 
of a = 12.5, subject to linearly x-direction polarized inci- 
dent microwave radiations of frequency a;/27r = 50 GHz 
having four different amplitudes Eis = 0.8,1.2,2.0 and 
3.5 V/cm at a lattice temperature of T = IK. For this 
GaAs system at 50 GHz aintor — 1-9 and ointra — 2.4. 
The intra-Landau level single-photon transitions are al- 
lowed when Wc/w > 2.4, yielding, at the high lUc/uj 
side, an absorption rate Sp somewhat larger, an electron 
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H =2500m /Vs a=12.5 
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E, = 0.1 V/cm 



FIG. 5: The magnetoresistivity Rxx vs ciJc/ii' for the same sys- 
tem as described in Fig. 3, subject to two incident microwave 
fields: frequency ti;/27r=30 GHz, amplitude Eis = 3 V/cm and 
frequency cj/27r=20 GHz, amplitude E^s ~ 2 V/cm. 



temperature Tg somewhat higher, and a SdHO damping 
stronger than those in the 100-GHz case (Fig. 1). On the 
other hand, at equivalent HF field strength the multi- 
photon processes are more important at lower frequency. 
This helps to enhance the absorption Sp in the range 
1.9 < Wc/w < 2.4, where the single-photon process is for- 
bidden and to increase the two-photon resonance in Sp 
and Tc around uj/uJc = 1-5,2.5 and 3.5 (see Sp and Tc 
curves corresponding to Eis ~ 3.5 V/cm in Fig. 3). The 
effect of the two-photon process can also be seen clearly 
in the Rxx-^^-^ l^c curves as shown in Fig. 4, where the 
Rxx curve of E-^s = 3.5 V/cm exhibits obvious shoulders 
around uj/uJc — 1.5,2.5 and 3.5, and the descends down 
around oj/uJc = 0.6. This kind of two-photon process was 
clearly seen in the experimentsi^ 

At even lower frequency, such as 30 GHz and 20 GHz, 
the ranges for intra-Landau level and inter-Landau level 
single-photon transitions overlap. The enhanced effect 
of the virtual photon process, together with enhanced 
multiphoton-assisted electron transition, pushes the re- 
sistivity Rxx remarkably down below the average of its 
oscillatory curve without radiation, resulting in a strong 
suppression of dissipative magnetoresistance across a 
wide magnetic field range as shown in Fig. 5, in agree- 
ment with experimental observations^i^ 



C. 150 and 280 GHz 

The radiation-induced SdHO modulation can be seen 
clearly in the low magnetic field region w/wc > 1 with 
higher radiation frequency. Figure 6 shows the calculated 
electron temperature Te and magnetoresistivity Rxx as 
functions of lo/oJc for a 2D system of electron density 
iVe = 3.0 X lO^^m-^, linear mobility /io = 2000mVVs 
and a = 3, subject to a 150-GHz microwave radiation of 




FIG. 6: The magnetoresistivity R^x and electron tempera- 
ture Te of a GaAs-based 2DEG with = 3.0 x lO^'^m-^, 
/io = 2000m^/Vs, and q = 3, subjected to 150-GHz linearly 
a;-polarized incident HF fields Eis smitut) of three different 
strengths Eis = 0.1,0.6 and 2 V/cm at T = 0.5 K. 



three different amplitudes E-^s = 0.1, 0.6 and 2 V/cm at a 
lattice temperature of T = 0.5 K. Low-power microwave 
illumination {E{s = 0.1 V/cm) already yields sufficient Te 
oscillation with maxima at lo/oJc = 1,2,3,4, giving rise 
to clear SdHO modulations having nodes at Te maxima. 
At higher microwave power {E^s = 0.6 V/cm) when the 
MIMO shows up, the Tc maxima gets higher, suppress- 
ing the SdHO in the vicinities of lo/lOc = 1,2,3,4, but a 
strong amplitude modulation of SdHOs is still seen. In 
the case of E^s = 2 V/cm, Rxx shows strong MIMO and 
the electron temperature further grows so that most of 
SdHOs almost disappear in the range of cj/oJc > 2. Note 
that the small Te peaks at w/wc = 1-5 and 2.5 are due to 
the absorption rate S'p maxima induced by two-photon 
processes, which gives rise to additional nodes in the Sd- 
HOs. 

Another example of the SdHO modulation appearing 
simultaneously with MIMO is plotted in Fig. 7, where the 
energy absorption rate S'p, the electron temperature Tc, 
and the magnetoresistivity Rxx are shown as functions of 
uj/uJc for a 2D system having an electron density of A^e = 
3.0 X 10^^ m~^, a linear mobility of ^jl^ = lOOOm^/Vs, and 
a broadening parameter of a = 2, subject to linearly x- 
direction polarized incident microwave radiations of fre- 
quency w/27r = 280 GHz and amplitude E^s = 3.5 V/cm. 
The energy absorption rate S'p has broad large peaks at 
uj/ujc = lj2,3,4, 5 (due to single-photon resonant pro- 
cess) and small peaks at uj/uJc — 1.5,2.5 (due to two- 
photon resonant process), giving rise to the oscillation 
of the electron temperature T^. One can clearly see the 
peaks of the electron temperature Te and the nodes of 
SdHO modulation at uj/oJc = 1,2,3,4 and 5, together 
with MIMOs. These are in agreement with the experi- 
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FIG. 7: The magnetoresistivity Rxx, electron temperature 
To and energy absorption rate Sp of a GaAs-based 2DEG 
with fj,o = lOOOm^/Vs and a — 2, subjected to a linearly 
a;-poIarized incident HF field of frequency 280 GHz and am- 
plitude Eis = 3.5 V/cm at T = 1 K. 

mental observation reported in Refini 

D. Discussion 

Note that in GaAs-based systems at a temperature 
around T ~ 1 K, LA phonons generally give larger con- 



tribution to the electron energy dissipation W than that 
from TA phonons and LO phonons are usually frozen. 
However, in the case of high radiation power or in the 
vicinity of a; ~ ujc, where the resonantly absorbed en- 
ergy can be relatively large and the electron temperature 
can rise up above 20 K, a weak emission of LO phonons 
takes place. Though at this temperature the number of 
excited LO phonons is still very small and their contri- 
bution to momentum relaxation (resitivity) is negligible 
in comparison with acoustic phonons, they can already 
provide an efhcient energy dissipation because each ex- 
cited LO phonon contributes a huge energy transfer of 
f^LO ~ 400 K. With a continuing rise of electron tem- 
perature the LO-phonon contribution increases rapidly. 
This effectively prevents the electron temperature from 
going much higher than 20 K, such that the Tc-vs-lUc/'-^ 
curve of large incident microwave power in Fig.l exhibits 
a flat top around lUc/^ = 1- 

In this paper, we did not consider the role of surface 
or interface phonons in the GaAs heterostructure. De- 
pending on sample geometry, the surface phonons may 
be important in dissipating electron energy thus decreas- 
ing the electron temperature. 
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